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The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are 
examined for the problem of hypersonic laminar shock/shock and shock/boundary layer 
interactions, where boundary layer separation is an important feature of the flow. Flow 
about two model configurations is considered, where both configurations (a biconic and a 
hollow cylinder-flare) have recent published experimental measurements. The computa- 
tions are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric 
time accurate code that incorporates many of the advances in DSMC over the past decade. 
The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test 
conditions with nitrogen as the test gas and the test models at zero incidence. Results 
presented highlight the sensitivity of the calculations to grid resolution, sensitivity to 
physical modeling parameters, and comparison with experimental measurements. Infor- 
mation is provided concerning the flow structure and surface results for the extent of 
separation, heating, pressure, and skin friction. 


Nomenclature 

a c,v surface accommodation for vibrational energy 
C f friction coefficient, TwRO.bpooV^) 

C h heat transfer coefficient, q w /(0-5pooV^o) 

C p pressure coefficient , (p w - Poo ) / (0 ■ 5poo ) 

d maximum double cone diameter, m 

L length of first cone projected on X-axis, m 

Moo free-stream Mach number 

m mass of a molecule, kg/m 3 

mcs mean collision separation distance, m 

mfp mean free path, m 

mpc molecules per cell 

n number density, m 3 

p pressure, N/m 2 

q wall heat transfer rate, kW/m 2 

R reattachment location, m 

Reoo unit Reynolds number, /.' N 1 \ /// x: . m _1 
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S separation location, m 

T temperature, K 

u„, velocity component normal to surface, m/s 

Up velocity component parallel to surface, m/s 
V s slip velocity at a surface, m/s 

Voo free-stream velocity, m/s 

x,y model coordinates, m 

p density, kg/m 3 

r viscous shear stress, N/m 2 

Subscripts 

R rotational temperature or reattacliment 
S separation 

T translational temperature 

V vibrational temperature 

W wall 

oo free stream 

Introduction 

H YPERSONIC flow interactions between shocks 
and separated flow regions produce very chal- 
lenging problems for both experimental measurements 
and numerical simulations. An understanding of such 
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interactions is critical to the design of hypersonic ve- 
hicles because these interactions can have a significant 
impact on both heating rates and vehicle aerodynam- 
ics, which in turn influence the requirements for vehicle 
control and thermal protection. Many experimental 
and theoretical studies dealing with shock/shock and 
shock/boundary layer interactions over compression 
ramps were conducted in the last half-century, and 
a recent review of some of the key findings for lam- 
inar flows is presented by Marini. 1 More recently, 
the NATO Research Technology Organization (RTO), 
under the coordination of Working Group 10, 2,3 has 
fostered a series of experimental and numerical studies 
focused on hypersonic flows. Two basic axisymmet- 
ric model configurations included in the RTO studies 
were a biconic and a hollow-cylinder flare. The current 
computational study focuses on the complex interac- 
tions resulting from hypersonic flows about these two 
configurations for experimental test conditions where 
separation and reattachment occur under laminar con- 
ditions. 

This paper presents computational results obtained 
with the direct simulation Monte Carlo (DSMC) 
method 4,5 for two of a larger set of test cases 6 that 
have been proposed for code validation studies. Both 
test cases are demanding for DSMC simulation since 
much of the flow domain is in the continuum regime; 
obvious exceptions are the flow near the sharp lead- 
ing edges and along the surface where velocity slip 
and temperature jump effects are significant. Further- 
more, both the locations and physical properties of the 
shocks and their interaction regions are very sensitive 
to the cell/grid resolution used in a simulation. The 
flow conditions simulated are those for two experimen- 
tal conditions performed by Holden and Wadhams: 6 
one for Mach 15.6 nitrogen flow about a 25°/55° bi- 
conic model in the Calspan-University at Buffalo Re- 
search Center (CUBRC) 48-Inch Shock Tunnel and the 
second from the measurements made for Mach 12.4 
nitrogen flow about a hollow cylinder-flare (30° flare) 
model in the CUBRC Large Energy National Shock 
(LENS) Tunnel. Details concerning the experimental 
measurements are given in Refs. 7 and 8, and a listing 
of the flow conditions for the two test cases considered 
in the present study is included in Ref. 9. 

Several DSMC studies have shown that the calcula- 
tion of such problems under continuum conditions can 
be very demanding in terms of computer resources, 
both in time and memory. In fact, results of the initial 
comparisons 10 of several DSMC calculations with mea- 
surements showed that the DSMC computations failed 
to predict the correct details of the separated region. 
Results of subsequent DSMC studies presented in Ref. 
9 suggest that the failure to predict the details of the 
separation region was due primarily to inadequate grid 
resolution. A second finding from the initial compar- 
isons of Harvey, et al. 10 was that both DSMC and CFD 


solutions differed substantially (as much as 33 percent, 
as discussed in Ref. 11) with the measured heating 
rates upstream of the separated region, that is, for lam- 
inar flow about both sharp 11 and blunted cones 12 at 
zero incidence. With differences of this magnitude, ad- 
ditional CFD studies 13 were conducted to re-examine 
the free-stream flow conditions for which the CUBRC 
measurements were made. The findings of the Can- 
dler, et al. 13 study for the CUBRC experiments showed 
that the free-stream conditions had a significant level 
of thermal nonequilibrium (initially assumed to be in 
equilibrium) resulting from vibrational freezing during 
the nozzle expansion, confirming the earlier assessment 
of Roy, et al. 12 The current calculations are made for 
one of the revised set of free-stream conditions 6 (hol- 
low cylinder flare, Run 11 conditions) and one for a 
biconic at a low density condition (Run 7). 

Experimental evidence 14 suggests that when there is 
significant vibrational nonequilibrium within the shock 
layer, as is the case in the current problems, the vibra- 
tional energy does not accommodate to the surface 
temperature when the molecules impact the surface. 
Calculated results for both test cases show a reduction 
of approximately 12 percent in surface heating rate 
when the surface is assumed to have zero rather than 
full thermal accommodation for the vibrational energy. 
Also, the calculations show that the assumption con- 
cerning surface vibrational energy accommodation has 
no impact on surface pressure or the extent of sepa- 
ration. Comparisons with the CUBRC measurements 
for heating and pressure show generally good agree- 
ment, with differences in heating on the order of 15 
percent or less, much less for these two test cases when 
full thermal accommodation is assumed. The rela- 
tionship between the slip velocity and the region of 
separation is investigated by calculating the sensitiv- 
ity of the slip velocity to numerical factors and specific 
physical factors associated with the gas-surface inter- 
actions. Results for this phase of the study indicate 
that the magnitude of the slip velocity and the ex- 
tent of separation can be significantly altered by the 
fraction of specular reflection associated with the gas- 
surface interactions, and the slip velocity cannot, by 
itself, be used to correlate the extent of separation. 

DSMC Program 

The DSMC program used in the current study is the 
DS2V program of Bird, 4,5 a general 2D/axisymmetric 
code that provides both time accurate unsteady flow 
and time-averaged steady flow simulations. Molec- 
ular collisions are simulated with the variable hard 
sphere molecular model. Energy exchange between 
kinetic and internal modes is controlled by the Larsen- 
Borgnakke statistical model. 15 For the present study, 
the simulations are performed by using nonreacting 
gas models while considering energy exchange between 
translational, rotational, and vibrational modes. The 
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nitrogen molecular gas constants used in the current 
study are those given in Ref. 4. Also, a rotational 
relaxation collision number of 5 and a temperature de- 
pendent vibrational collision number (Eq. 6.53 of Ref. 
4) were used. 

The DS2V program has the option of specifying 
multiple free-stream temperatures, as is appropriate 
for the test cases currently investigated. The model 
surface is assumed to have a specified constant tem- 
perature, and diffuse reflection is assumed for the 
gas-surface interactions unless noted otherwise. Two 
options are exercised concerning surface thermal ac- 
commodation: one in which all energy modes are fully 
accommodated and one in which all energy modes, ex- 
cept for vibration, are fully accommodated (a c,v = 0). 

The DS2V program has taken advantage of several 
important developments 5 in DSMC methodology over 
the past decade, two of which will be briefly described. 
One is the use of transient subcells, in which a tran- 
sient background grid is applied one cell at a time 
within the collision routine to promote nearest neigh- 
bor collisions. The resolution of the transient back- 
ground grid is a function of the number of molecules in 
the cell, that is, approximately one simulated molecule 
within each grid subcell. The ratio of the mean colli- 
sion separation between collision partners to the local 
mean free path ( mcs/mfp ) should be less than one, and 
this parameter is currently being used as an indicator 
of the resolution achieved in a given simulation. One 
option for achieving a small value for this parameter is 
to increase the number of simulated molecules, which 
can be exercised any time during the simulation. Such 
an approach is demonstrated in the current study. 

A second new feature of the current program is the 
manner in which the time step is established for un- 
coupling the molecular motion and collisions. Rather 
than use a single time step over the whole flow field 
or multiple subdomains, with each domain having a 
different time step, new procedures are used in which 
the time step varies with every molecule and every cell. 
These procedures, described in Ref. 5, keep all the time 
variables concurrent with the overall flow time. The 
overall flow time step is advanced in very small time 
increments (a fraction of the local mean collision time, 
which is kept for each cell). For an overall time step, 
only a small fraction of the molecules is moved, and 
only a small fraction of the molecules experiences col- 
lisions. 

The flow-field cells are created from rectangular 
structures that consist initially of divisions, but from 
elements after the first grid adaption. After adap- 
tion, the cells are not necessarily rectangular because 
they are formed from one or more elements. The 
computational domain which is initially divided into 
a Cartesian network of cells with a specified num- 
ber of divisions in each of the coordinate directions 
is bounded by user specified limits. If one is simulat- 


ing a steady flow problem, the usual procedure is to 
advance the simulation until the flow achieves or ap- 
proaches a steady state and then do a grid adaption. 
Adaption is achieved by creating new cells made up of 
one or more Cartesian elements (each division contains 
a fixed number of elements) so that each cell contains 
a specified number of simulated molecules. The tran- 
sient subcells are set for each cell whenever collisions 
are to be calculated and are then discarded. The use 
of the transient subcells leads to nearest neighbor col- 
lisions. 

Note that during a calculation interruption to ei- 
ther adapt the grid or change the number of simulated 
molecules, information at the time of interruption is 
retained, and the calculation continues to advance in 
time, with some time interval required before the per- 
turbation effects, caused by an interruption, have been 
assimilated. The cell structure is set largely by the 
DS2V program logic and relatively little input is re- 
quired of the user, which leads to a more user friendly 
program. Also, the new procedures are designed to 
be computationally fast, and the built-in flow visual 
displays provide instant access to the simulation pa- 
rameters and results as they evolve with time. 

An option that DS2V provides was exercised in the 
current study in which the surfaces may be set to move 
in the flow plane. Also, surfaces can be set to move 
normal to the flow plane, which permits the study of 
rotating flows for axisymmetric bodies. Cells for sam- 
pling properties along solid surfaces may be specified 
at any desired spacing. 

Results and Discussion 

Table 1 presents a summary of the free-stream con- 
ditions for the two test cases investigated. The current 
paper presents results calculated with the DS2V pro- 
gram for the biconic at Run 7 test conditions, as 
generated in the 48-Inch Shock Tunnel, and for the 
hollow cylinder-flare model at Run 11 test conditions, 
as generated in the LENS I Tunnel. The Mach and 
Reynolds numbers listed are those provided for the re- 
spective test facilities and are based on the free-stream 
translational temperatures. The free-stream mean free 
path is approximately 0.219 mm and 0.084 mm for Run 
7 and Run 11 conditions, respectively. 

Biconic Calculations for CUBRC Run 7 Test 
Conditions 

The geometry and dimensions for the sharp 25°/55° 
double cone model used in the current calculations are 
shown in Fig. 1. Reference 8 provides information 
regarding the experimental procedures and instrumen- 
tation that was used to obtain the heat transfer and 
pressure measurements that are compared with the 
present computations. 

Table 2 presents a summary of some of the numer- 
ical parameters used in the calculations, along with 
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Table 1 Free-stream and surface conditions for biconic and hollow cylinder-flare calculations. 


Condition 

Voo , rn/s 

?ioo, m 3 

Too,T, K 

Poo, kg/m 3 

Poo, N/m 2 * 

Mqo 

Re oo, m 1 

Gas 

T W , K 

48-Inch Run 7 

2 073 

3.779 x 10 21 

42.6“ 

1.757 x 10 -4 

2.23 

15.6 

137467 

Ni 

297.2 

LENS Run 11 

2 484 

1.197 x 10 22 

95.6 6 

5.566 x 10 -4 

15.86 

12.4 

208333 

No 

297.2 


" Too ,T = Toc.K = 42.6 K and T*y = 1986 K. 
b ^oc,T = TcK'R = 95.6 K and T MiV = 2487 K. 


Table 2 DSMC numerical parameters 0 and results for sharp double cone (25°/55°), Run 7 test conditions. 


Solution 

A time, ms 

mpc 

Molecules 

mcs /mfp 

Vib. Accom. 

xs, mm 

xr : mm 

Ax, mm 

aO 

1.27 to 1.91 

10 

0.514E+6 

0.944 

0.0 

82.8 

99.4 

16.6 

bO 

3.30 to 3.65 

10 

1.980E+6 

0.441 

0.0 

81.2 

100.4 

19.2 

cO 

3.84 to 3.97 

40 

8.136E+6 

0.214 

0.0 

80.3 

100.7 

20.4 

dO 

3.97 to 4.41 

40 

8.114E+6 

0.213 

0.0 

80.4 

100.8 

20.4 

al 

1.03 to 2.02 

10 

0.528E+6 

0.946 

1.0 

82.6 

99.6 

17.0 

bl 

3.30 to 3.65 

10 

1.980E+6 

0.441 

1.0 

81.1 

100.4 

19.3 

cl 

3.83 to 3.96 

40 

8.133E+6 

0.214 

1.0 

80.4 

100.8 

20.4 

dl 

3.96 to 4.40 

40 

8.114E+6 

0.213 

1.0 

80.4 

100.8 

20.4 

el 

4.81 to 4.93 

40 

16.482E+6 

0.147 

1.0 

80.1 

101.0 

20.9 


c Divisions = 151 x 181, elements per division = 20 x 20, and radial weighting factor = 100. 


193.7 



Fig. 1 CUBRC sharp double cone model where x 
is measured from the vertex. 


results for the locations and extent of the separation 
region. Surface sampling cells were specified as fol- 
lows: 100 along the fore-cone, 150 along the second 
cone, and 15 along the model extension. Figures 2 
through 15 present details concerning the calculated 
flow structure and surface results. Two separate se- 
ries of calculations are made where the surface thermal 



Fig. 2 Density contours from solution el for 
CUBRC Run 7 conditions. 

accommodation for the vibrational energy is the only 
variable, zero accommodation for the lettered solutions 
with a “0” extension, and full accommodation for the 
lettered solutions with a “1” extension, as listed in 
Table 2. Results for these time accurate solutions are 
included for flow times of approximately 1 to 5 ms. 

Figures 2 through 6 present flow-field contours that 
describe many of the flow features for the current prob- 
lem and the computational domain (Fig. 2) used to 
perform the numerical simulations. Figure 7 presents 
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E 


0.065 

0.060 

0.055 

0.050 

0.045 

0.040 

0.035 



x, m 


Fig. 3 Density contours in separation region for Fig. 5 Overall kinetic temperature contours and 
Run 7 conditions, solution el. streamlines in separation region for Run 7 condi- 

tions, solution el. 



Fig. 4 Pressure contours and streamlines in sepa- 
ration region for Run 7 conditions, solution el. 

contours for the ratio of the mean collision separation 
distance to the mean free path ( mcs/mfp ), a param- 
eter used to measure the quality of the simulation 
where it is desirable that the parameter be less than 
one throughout the computational domain. Further- 
more, Figs. 8 through 15 provide calculated surface 
distributions for velocity slip, temperature jump, heat- 
ing, pressure, and friction, along with some insight 
on how these quantities are influenced by computa- 
tional resolution, how the calculated data compare 
with the CUBRC surface heating and pressure mea- 
surements, and how the current computations com- 
pare with those from another DSMC program called 
SMILE 16 (Statistical Modeling In Low-density Envi- 
ronment) . 

Closer examination of the data presented in Figs. 2 


E 


0.065 

0.060 

0.055 

0.050 

0.045 

0.040 

0.035 



0.08 0.09 0.10 0.11 


x, m 


Fig. 6 Vibrational temperature contours and 
streamlines in separation region for Run 7 condi- 
tions, solution el. 

through 15 reveals several points of interest. First, the 
application of the DS2V code provides results that are 
in good (less than 10 percent) to fair (10 to 15 per- 
cent difference) overall agreement (Figs. 10 through 
15) with the experimental data. Second, the thermal 
nonequilibrium free-stream condition persists within 
the shock layer (compare Figs. 5 and 6). Only within 
the boundary layer and near the surface is there evi- 
dence (Fig. 6) of significant depletion or equilibration 
of the vibrational energy, and only for the solution with 
full thermal accommodation at the surface. When the 
calculation is performed with no surface accommoda- 
tion of the vibrational energy (solution dO), the vibra- 
tional temperature remains above 1925 I< throughout 
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Fig. 7 Ratio of mean collision separation to mean 
free path (mcs/mfp) contours and streamlines in 
separation region for Run 7 conditions, solution el. 



Fig. 9 Calculated surface distributions for heating, 
pressure, and friction for Run 7 conditions, solution 
el. 



Fig. 8 Calculated surface distributions for velocity 
slip and temperature jump, both translational and 
rotational, for Run 7 conditions, solution el. 

the flow field, essentially frozen at the elevated free- 
stream condition. Calculations performed with the 
same free-stream conditions, except that the vibra- 
tional temperature was set to that of the other temper- 
ature modes (42.6 K), resulted in lower heating rates 
(not shown) but showed no perceptible differences for 
surface pressure or friction. Third, the current time ac- 
curate simulation suggests that this laminar flow test 
case, with boundary layer separation, stabilizes rather 
quickly, on the order of a millisecond. The longer run 
time used in the present study occurred because sev- 
eral numerical perturbations (grid adaption or increase 
in the number of simulated molecules) were imposed 
as the simulation evolved in time. The simulations 


Fig. 10 Calculated and measured heating rates for 
Run 7 conditions — effect of computational resolu- 
tion. 

were monitored to ensure that a steady state was 
achieved, followed by the generation of a significant 
time-averaged sample for the results listed in Table 2. 
The steady state aspect of the simulations is demon- 
strated in Table 2 for the cO/1 and dO/1 solutions, 
both with and without surface vibrational accommo- 
dation, where sequential time-averaged solutions yield 
essentially the same results for the separation region. 
Fourth, the current simulation approach has used com- 
putational cells that are coarse in relation to the local 
mean free path, but the simulation has used the tran- 
sient subcell feature of the DS2V code, along with a 
large number of simulated molecules per cell, to pro- 
mote small values of the local mcs/mfp parameter. For 
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Vib. Accom. = 1 .0, solution dl 
Vib. Accom. = 0.0, solution dO 




Fig. 11 Calculated and measured surface pressures 
for Run 7 conditions — effect of computational res- 
olution. 


Fig. 13 Calculated and measured surface pressures 
for Run 7 conditions — effect of vibrational surface 
accommodation. 


DS2V, Surf, vibrational Accom. = 0.0, solution dO 
SMILE, No vibrational energy 




Fig. 12 Calculated and measured heating rates 
for Run 7 conditions — effect of vibrational surface 
accommodation. 

the most resolved simulation, solution el, the mean 
value of mcs/mfp was 0.15, with local values less than 
0.37 at the location of separation and with no value 
exceeding 0.72. 

Additional numerical studies would be useful in es- 
tablishing bounds on the computational merit param- 
eter mcs/mfp to ensure good engineering results for 
general applications. Obviously, if the parameter is 
sufficiently small, the results should be accurate but 
potentially demanding computationally. Also, a small 
mean value of the merit parameter is at best only a 
necessary condition for a good simulation, but not a 
sufficient condition because of potentially large gradi- 
ents, whereas local values of the merit parameter can 


Fig. 14 Calculated and measured heating rates for 
Run 7 conditions — results from two DSMC codes. 

be much larger than the mean value. 

The biconic test cases produce large flow-field gra- 
dients due to shock interactions and substantial am- 
plifications of surface loads caused by shock/boundary 
layer interactions, as is clearly evident for the CUBRC 
Run 7 test case. Densities that equal 30 times the free- 
stream value are evident within the shock layer (Fig. 
3), and even larger densities (50 times free stream) 
occur along the second cone surface. The peak pres- 
sure amplification occurs off the surface (Fig. 4) with a 
magnitude that equals 670 times the free-stream value, 
while the maximum overall kinetic temperature (Fig. 
5) is 2050 K. 

In Fig. 8, the calculated surface distributions for ve- 
locity slip and temperature jump, both translational 
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DS2V, Surf, vibrational Accom. = 0.0, solution dO 


SMILE, No vibrational energy 



Fig. 15 Calculated and measured surface pressures 
for Run 7 conditions — results from two DSMC 
codes. 

and rotational, are presented. The maximum values 
for these quantities occur at or near the cone apex, 
then decrease with increasing distance downstream, 
decrease significantly in the separation region with a 
reversal in the sign of the slip velocity, and increase 
with the flow expansion at the aft end of the second 
cone. Similar results to those shown in Fig. 8 will be 
discussed in more detail in the next section. 

The calculated surface distributions for heating rate, 
pressure, and friction are presented in Figs. 9 through 
15, in which Fig. 9 presents a summary of heating, 
pressure, and friction results (in coefficient form) for 
the most resolved calculation (solution el), 16.5 mil- 
lion simulated molecules with a mean mcs/mfp value 
equal to 0.147. The maximum local values for the 
mcs/mfp parameter never exceeded 0.72, which oc- 
curred along the second cone. Figures 10 and 11 show 
the effect of reducing the mean value of mcs/mfp (first 
by a factor of 2 [increasing the number of molecules by 
a factor of 4] and then by a factor of 1.41) on heating 
and pressure distributions, respectively. The results 
included in Table 2 indicate the corresponding effect 
of reducing the magnitude of the mcs/mfp parameter 
from about 0.95 to 0.15 on the locations and extent 
of separation, with the extent of separation increas- 
ing with increased resolution. As shown in Figs. 10 
and 11, the agreement between the finest resolution 
simulation (with a c,v = 1-0) and the measurements 
is reasonably good, with the calculation predicting a 
slightly delayed separation relative to the experimental 
value. 

The impact of the surface boundary condition for 
the vibrational energy accommodation is clearly evi- 
dent on the surface heating, as is demonstrated in Fig. 
12 where the heating rate is reduced by approximately 
12 percent for zero accommodation, as compared to 


full accommodation. The impact of the vibrational 
energy accommodation boundary condition on surface 
pressure (Fig. 13) and the locations and extent of 
separation (Table 2) are negligible. The comparisons 
shown in Figs. 12 and 13 are made for the same grid 
resolution, solutions dl and dO. The calculated heat- 
ing rates are in better agreement with the CUBRC 
measurements for the assumption of full thermal ac- 
commodation and the same is shown later for the Run 
11 test case; therefore, one might be inclined to say 
that the full accommodation of the vibrational energy 
mode is more appropriate. However, results from addi- 
tional calculations for only the forecone of the biconic 
model at other CUBRC test conditions (Runs 28 and 
35 as discussed in Ref. 17) show a different trend where 
the zero accommodation calculation for the vibrational 
energy mode provides better agreement with measure- 
ments. Therefore, it appears that the variability in the 
measured heating rate data, with respect to the DS2V 
calculations, is such that one is not able to resolve 
a 12 percent effect in heating; that is, to determine 
if either of the vibrational accommodation boundary 
conditions is appropriate. 

The final comparison for the Run 7 calculations are 
the current DS2V results with those Markelov ob- 
tained by using the SMILE DSMC code, the same 
results that were included in a recent review paper 
by Harvey (Ref. 9). The SMILE calculations were 
made without accounting for the vibrational energy 
(vibration turned off), which should be very similar 
to the DS2V results with vibration included and zero 
accommodation for the vibrational energy at the sur- 
face. Thus the SMILE results are compared with data 
from solution dO. The two calculations for surface 
heating and pressure distributions are in very close 
agreement, as shown in Figs. 14 and 15, with the 
SMILE results showing slightly better agreement with 
the experimental measurements, particularly for the 
extent of separation ( xs = 79.6 mm, xr = 101.1 mm, 
and Ax = 21.5 mm). In the SMILE calculation, a par- 
allel code implementation with 80 million simulated 
molecules was used, approximately an order of mag- 
nitude grater than that used in the single processor 
DS2V calculation, solution dO. 

Relationship between the Slip Velocity and the 
Region of Separation 

The slip velocity V s at a surface element may be 
determined from the incident and reflected velocities of 
the molecules that strike the element. If m is the mass 
of a molecule, u p is the velocity component parallel to 
the surface, and u n is the velocity component normal 
to the surface, 

v s = Yl ( mu p / 1 Un i) / ( m / 1 Un i) w 

Note that the sums are taken over both the incident 
and the reflected molecules. The expression for the slip 
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velocity given by Eq. 1 differs from the standard defi- 
nition of velocity in a spatial or volumetric element by 
the addition of the terms containing the normal veloc- 
ity component. These additional terms transform the 
fluxal quantities in the surface samples to the spatial 
quantities in the gas immediately above the surface. 
Similar transformations have been made to the spatial 
definitions of translational and rotational temperature 
in order to determine the temperature jumps. 

The slip velocity in a simulation is affected by the 
DSMC computational merit parameter mcs/mfp , as 
well as by physical factors such as the gas-surface in- 
teraction through the inclusion of a fraction of specular 
reflection at the surface. It may also be affected by 
moving the surface in its own plane in the stream-wise 
direction. To highlight the effect of the slip velocity 
on the extent of separation, the DS2V program was 
used to make a set of runs (Tables 3 and 4) for the 
Run 7 conditions that are outlined in Table 1. Each 
calculation employed in this phase of the study had ap- 
proximately 2.5 million simulated molecules, and the 
final flow sample was from 1.5 to 2 ms. A negative 
surface velocity is towards the apex of the cone, and 
the molecules are reflected from the diffusely reflecting 
surface then acquire a velocity increment that reduces 
the slip velocity. 

The slip velocity is a maximum at the apex of the 
cone and declines with distance from the apex (similar 
to the results shown in Fig. 8 for solution el, Table 
2). Table 3 shows the slip velocity at a location just 
upstream of the earliest separation point, and there is 
an almost linear relationship between the slip velocity 
and the x coordinate of the separation point. The 
separation point moves upstream by 0.2 mm for each 
1 m/s reduction in the slip velocity. Similarly, the 
reattachment point moves downstream by 0.075 mm 
for each 1 mm reduction in the slip velocity at x = 60 
mm. 

The 2.5 million simulated molecules for the slip 
sensitivity calculations lead to an mcs/mfp ratio of 
0.6 in the vicinity of the separation point, where the 
corresponding value was 0.23 for the most resolved 
calculation listed in Table 2 (solution el, with 16.5 
million simulated molecules). A calculation with ap- 
proximately 25 million molecules would be required to 
bring the mcs/mfp ratio down to the recommended 
maximum value of 0.2, since the mcs/mfp ratio de- 
creases by a factor of two for a factor of four increase 
in the number of simulated molecules. While this has 
not been possible in the present study, calculations 
of this quality have been made of the flow over the 
fore-cone. The slip velocity at x = 60 mm in these 
calculations was 51 m/s and, assuming that the slip ve- 
locity in an imperfect calculation (too large a mcs/mfp 
ratio, which produces too larger a slip velocity) has the 
same effect on the separation as the in-plane surface 
velocity, the converged separation point would be at x 


= 79.5 mm and the reattachment point at x = 100.9 
mm. These values are in good agreement with both 
the measured values and the nearly grid-converged cal- 
culations presented in the previous section. 

The preceding comparisons and discussion indicate 
that changes in the slip velocity have the same effect on 
the separation irrespective of whether the changes are 
produced by the calculation’s spatial resolution or by 
in-plane surface movement. It is tempting to extend 
this finding to other factors that affect the slip velocity, 
but calculations that employ a combination of diffuse 
and specular reflection at the surface have shown that 
such an approach would be invalid. 

A small fraction of specular reflection has little effect 
on the extent of separation but causes a large increase 
in the velocity slip at the surface, as indicated in Table 
4. If there was a positive in-plane surface velocity suf- 
ficient to increase the slip velocity at x = 60 mm to 135 
ms/, as in the 0.33 specular case, the separation would 
be almost completely suppressed, as demonstrated in 
Table 3. It is clear that the slip velocity cannot, by 
itself, predict the extent of separation. 

The Mach number at the surface in the 0.66 specular 
case is supersonic and, if the flow was in equilibrium, 
there could be no upstream propagation of distur- 
bances and no separation. However, the velocity dis- 
tribution at the surface is a near bimodal combination 
of the distributions that correspond to the diffusely re- 
flected molecules and the stream molecules. It is the 
diffuse-like component of the velocity distribution that 
permits the upstream propagation of disturbances. 

Traditional discussions of rarefied gas flows assume 
the existence of a “slip flow regime” such that the 
flow may be validly described by the Xavier Stokes 
equations with slip boundary conditions at surfaces. 
On the other hand, it has been pointed out that the 
Navier-Stokes equations do not contain the physical 
phenomena that produce the velocity slip and temper- 
ature jump. The biconic flows with separation that 
are studied in this paper demonstrate that these mis- 
givings are justified. Depending on factors that do not 
enter the slip model, very different flows may be as- 
sociated with a given level of velocity slip. The “slip 
flow regime” is indeed a myth! 

A small fraction of specular reflection also has little 
effect on the net heat transfer rate for the Run 7 con- 
ditions. As the specular component of reflection was 
increased from 0 to 100 percent for the cases listed in 
Table 4, the corresponding heat transfer rates at the 
60 mm location were 33, 32, 29, and 0, respectively. 

Hollow Cylinder-Flare Calculations for CUBRC 
Run 11 Test Conditions 

The CUBRC Run 11 test case for the hollow 
cylinder-flare model was included in the present study 
because it has been investigated extensively with both 
CFD and DSMC methods for the initially reported 9 10 
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Table 3 Effect of an in-plane stream-wise surface velocity on the extent of separation. 


Surface Speed, m/s 

xsi mm 

xr, mm 

A a;, mm 

Y s at x = 60 mm, m/s 

-100 

68.0 

106.0 

38.0 

5 

-50 

76.3 

102.8 

26.5 

32 

-25 

78.5 

101.4 

22.9 

46 

0 

80.9 

100.4 

19.5 

58 

25 

82.8 

98.5 

15.7 

72 

50 

84.7 

97.0 

12.3 

86 

Table 4 Effect of gas- 

-surface interaction 

model on 

the extent of separation. 


Fraction of Specular Reflection 

xs, mm 

xr, mm 

Ax, mm 

V s at x = 60 mm, m/s 

0 

80.9 

100.4 

19.5 

58 

0.33 

81.8 

100.3 

18.5 

135 

0.66 

88.2 

95.3 

7.0 

717 

1 

- 

- 

0.0 

1840 


test conditions, has a small but distinct separation 
region, and has been examined recently with DSMC 
computations 9,18 as the free-stream conditions were 
redefined 6 to account for vibrational nonequilibrium 
effects. Figure 16 presents details of the model config- 
uration used in the present calculations. The hollow 
cylinder has a sharp leading edge with a bevel angle of 
10° and is aligned parallel to the oncoming flow. The 
compression flare is inclined 30° to the cylinder. Sur- 
face sampling cells were specified as follows: 600 along 
the flare, 250 along the top surface of the cylinder, and 
20 along the bevel leading edge. 

An abbreviated presentation of the calculated re- 
sults for the Run 11 test case is given in Table 5 
and Figs. 17 through 19, highlighting the comparison 
with the experimental measurements for surface heat- 
ing and pressure. Figure 17 shows the computational 
boundary used in the present calculations, which in- 
cludes a portion of the beveled leading edge. Also 
shown in Fig. 17 are selected density contours nor- 
malized by the free-stream value, where the maximum 
value is 26. 

Similar to the Run 7 solutions, two separate sets of 
calculations were made for Run 11 conditions where 
the only variable was the surface accommodation for 
the vibrational energy. As listed in Table 5, the set 
of solutions with zero surface accommodation for the 
vibrational energy is denoted by a letter followed by 
“0” and the solutions with full accommodation are de- 
noted by a letter followed by “1”. For each set of 
calculations, some of the numerical parameters and 
results for the locations and extent of separation are 
included in Table 5 as the solutions were advanced in 
time with various perturbations accounting for grid 
adaption and increases in the number of simulated 
molecules. The steady state aspect of the simulations 


is demonstrated in Table 5 for the g0/l and h0/l so- 
lutions, both with and without surface vibrational ac- 
commodation, where sequential time-averaged results 
provide essentially the same results for the locations 
and extent of the separation region. The best spatially 
resolved solutions (jO and jl) had global mean values 
for mcs/mfp of 0.31 and maximum values of 0.25 and 
0.75 at separation and reattachment, respectively. The 
maximum value of mcs/mfp was 1.2, and this value oc- 
curred along the surface in the region downstream of 
reattachment where the shock interactions produce a 
very thin shock layer. 

Results calculated for the surface heating rate and 
pressure distributions, along with the location and ex- 
tent of separation, are included in Fig. 18. Agreement 
between the calculated and the measured results are 
generally good for both the heating rate and pressure 
distributions (solution jl, Table 5), particularly the 
heating rate. For this test case, the calculated location 
of separation is slightly upstream of the results inferred 
from the heating and pressure measurements. Also ev- 
ident in Fig. 18 are the rarefaction effects at the hollow 
cylinder leading edge, where both the heating rate and 
pressure increase as the leading edge is approached, 
reach their respective maximum and then decrease to- 
ward their free molecular values, failing to achieve the 
free molecular values because of upstream influence 
from the leading edge. Along the cylinder, the veloc- 
ity slip and temperature jump are substantially larger 
than those experienced by the biconic model for Run 
7 conditions (Fig. 8), with leading edge values of ap- 
proximately 875 m/s for the velocity slip and 1400 K 
for the translational temperature jump. 

The calculated effect of surface thermal accommo- 
dation is addressed in the results presented in Fig. 19. 
The heating rate distribution, assuming full thermal 
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accommodation (blue curve, solution jl in Table 5), 
is about 12 percent higher than that obtained for the 
assumption that the surface vibrational energy is not 
accommodated (red curve, solution jO in Table 5). As 
discussed in Ref. 14, there are experimental results for 
vibrationally nonequilibrium flows which suggest that 
the surface vibrational energy accommodation is very 
small and that a zero accommodation should be a real- 
istic boundary condition for the present problems. As- 
suming full accommodation, the current calculations 
provide good agreement with the heating measure- 
ments, but only fair agreement with the assumption 
of zero accommodation for the vibrational energy at 
the surface. As was the case for the Run 7 calcula- 
tions, the effect of the surface boundary condition on 
the separation region and the pressure distribution was 
negligible. 



Fig. 18 Calculated and measured surface heating 
rate and pressure distributions for Run 11 condi- 
tions, solution jl. 


Fig. 16 CUBRC hollow cylinder-flare model (di- 
mensions in mm) where x is measured from leading 
edge. 



Fig. 17 Computational domain and density con- 
tours for Run 11 conditions, solution jl. 



Fig. 19 Calculated effect of vibrational energy 
surface thermal accommodation on heating rate 
distribution for Run 11 conditions, solutions jO and 

jl- 

Concluding Remarks 

Results of a computational study are presented for 
hypersonic flows about biconic and hollow cylinder- 
flare models at zero incidence. The free-stream con- 
ditions are those given for experiments conducted by 
CUBRC at Mach 15.6 and 12.4 in nitrogen in two 
shock tunnels, the 48-Inch Tunnel and the LENS I 
Tunnel, respectively. Computations are made with the 
DSMC method by using a relatively new and evolv- 
ing time-accurate program under development by Bird 
(the second author), called DS2V. The focus of the 
current study is to examine the effectiveness of the 
DS2V code in predicting complex flows that result 
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Table 5 DSMC numerical parameters^ and results for a hollow cylinder-flare, Run 11 test conditions. 


Solution 

A time, ms 

mpc 

Molecules 

mcs / mfp 

Vib. Accom. 

xs, mm 

xr : mm 

Arc, mm 

fO 

0.92 to 1.02 

10 

1.707E+6 

0.781 

0.0 

93.6 

122.1 

28.5 

gO 

1.24 to 1.57 

10 

1.705E+6 

0.782 

0.0 

91.1 

122.8 

31.7 

hO 

3.64 to 4.29 

10 

1.639E+6 

0.768 

0.0 

90.8 

122.5 

31.7 

iO 

4.49 to 4.79 

31 

5.200E+6 

0.440 

0.0 

88.8 

124.9 

36.1 

jo 

4.91 to 5.11 

40 

10.339E+6 

0.306 

0.0 

86.1 

126.3 

40.2 

fl 

0.90 to 1.01 

10 

1.682E+6 

0.775 

1.0 

94.2 

121.2 

27.0 

gl 

1.14 to 1.92 

10 

1.706E+6 

0.785 

1.0 

90.8 

122.2 

31.4 

hi 

3.02 to 3.67 

10 

1.641E+6 

0.770 

1.0 

90.4 

122.7 

32.3 

il 

3.90 to 4.12 

31 

5.106E+6 

0.436 

1.0 

88.6 

124.9 

36.3 

jl 

4.71 to 4.85 

40 

10.337E+6 

0.307 

1.0 

86.2 

126.3 

40.1 


d Divisions = 223 x 91, elements per division = 20 x 20, and no radial weighting factor. 


from shock/shock and shock/boundary layer interac- 
tions and problems where the calculated results are ex- 
tremely sensitive to the numerical parameters used in 
the simulation. The overall comparisons of the present 
computations with both experimental and other com- 
putational results are shown to be good. 

The results obtained for the CUBRC Run 7 test 
case (biconic model) suggest the following: The cur- 
rent time-accurate calculations, made with a coarse 
cell structure but with a large number of molecules per 
cell (fairly small values of the mean collision separa- 
tion to mean free path ratio (mcs/mfp)), give generally 
good agreement with the CUBRC data for surface 
heating and pressure. In particular, the heating rate 
comparisons show good agreement for the most refined 
calculations, with the assumption of full surface ac- 
commodation for the vibrational energy, but only fair 
agreement for the assumption of no surface accommo- 
dation of the vibrational energy, which produces a 12 
percent reduction in the heating relative to the full ac- 
commodation assumption. Comparison of the present 
DS2V results with those obtained with the SMILE 
DSMC program shows very good agreement. Also, the 
thermal nonequilibrium free-stream condition, with a 
highly elevated free-stream vibrational temperature, 
persists throughout the shock layer for the assump- 
tion of zero surface accommodation of the vibrational 
energy. When a full thermal accommodation surface 
boundary condition is imposed, the vibrational energy 
within the shock layer experiences a significant deple- 
tion or adjustment adjacent to the surface, producing 
approximately a 12 percent increase in the surface 
heating rates. 

For the CUBRC Run 11 test case (hollow cylinder- 
flare model), the computed and measured surface re- 
sults for heating rate and pressure showed the same 
general trends evident for the Run 7 case; good agree- 
ment for heating rates with full surface accommoda- 
tion but only fair agreement with zero surface accom- 
modation for the vibrational energy (also, a 12 percent 


reduction in surface heating rates). 

Even though the present comparisons of measured 
and calculated heating rates suggest that the full sur- 
face accommodation is the appropriate boundary con- 
dition, other calculations presented in Ref. 17 using the 
same DS2V program for only the fore-cone of the bi- 
conic model and at higher density conditions show the 
opposite trend; that is, the use of the zero accommoda- 
tion for the vibrational energy provides better agree- 
ment with the CUBRC measurements. Consequently, 
the variability evident in the heating rate measure- 
ments, when compared to the present calculations, are 
such that it is difficult to isolate a 12 percent effect on 
heating rates based on the comparisons for four dif- 
ferent test cases; that is, which boundary condition 
is more appropriate for the surface accommodation of 
the vibrational energy. 

For the Run 7 test conditions, calculations are made 
that examine the relationship between the slip velocity 
and the region of separation for a diffuse surface, where 
it is shown that a correlation can be established for 
relating the size of the separation region with the slip 
velocity, but such a correlation becomes invalid when 
the physics of the gas-surface interactions includes 
some small fraction of specular reflection. Therefore, 
depending on the physics included in a slip model, very 
different flows may be associated with a given level of 
velocity slip. 
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